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Asymptotic behavior of solutions to space-time fractional 

diffusion equations 
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Abstract 

This article discusses the analyticity and the long-time asymptotic behavior of solutions 
to space-time fractional diffusion equations in R d . By a Laplace transform argument, we 
prove that the decay rate of the solution as t —> oo is dominated by the order of the 
time-fractional derivative. We consider the decay rate also in a bounded domain. 
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1 Introduction and main results 

Letting 0 < a < 1 and 0 < 0 < 2 be fixed constants, we first deal with the Cauchy problem for 
the space-time fractional diffusion equation 

| d?u(x,t) = -(-A )%u(x,t)+p(x)u{x,t), ^ 

\ u(x, 0) = a(x), x £ t > 0. 

Here we point out that the space-time fractional diffusion equation under consideration is 
obtained from the standard diffusion equation by replacing the second-order space derivative 
with a Riesz derivative of order /? £ (0,2] defined by 

(-A )Su(x,t):= f |^u(£,f)e 2 ^df, 

JR d 

where u denotes the Fourier transform u(£,t) := f Rd u{x,t)e~ 2m ^' x dx, and the first-order time 
derivative with a Caputo derivative of order a £ (0,1) that is usually defined by the formula 

d?u(x,t) := —2- t( (t-T)~ a ^(x,T) dr, t > 0. 

T(1 - a) J o dr 

From the physical view point of the time-fractional and space-fractional derivatives respectively, 
we see that replacing a time derivative with its fractional counterpart basically introduces 
memory effects into the system and then makes the process non-Markovian, while a replacement 
of a space derivative introduces non-local effects. The fundamental solution (for the Cauchy 
problem) of space-time fractional diffusion equation can be interpreted as a probability density 
evolving in time of a peculiar self-similar stochastic process. See, e.g., [5], [16] and the references 
therein. 
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In some recent publications such as e.g. [2] and [7], the space-time fractional diffusion 
equations on bounded domain are investigated. The space-time fractional diffusion equations 
are usually defined as 

dfu =—A^u+p(x)u, is!l, t > 0, ( 1 . 2 ) 

where 0<a<l,0</3<2, flcM d isa bounded domain with smooth boundary dtt and 
A := —A with D(A) = H 2 (Tl) n Hq(Q). See e.g. [18j for the definition of A%. 

For a = 1 and 0 = 2, the equations ED and ED turn to be the classical parabolic 
equations. For such kind of equations we refer Friedman [3] or Pazy |18| for the long-time 
asymptotic estimate of the solutions. We also refer Yamamoto [27] where several kinds of 
Carleman estimates are established and applied to some inverse problems. 

For a = 1 and 0 G (0,2), it becomes the space fractional diffusion equation. See, e.g., [1], 
m and the references therein. 

For a € (0,1) and 0 = 2, it is the so-called time-fractional diffusion equation. For such a 
kind of equations there are a large and rapidly growing number of publications. In the 
wcll-posedness and the long-time asymptotic behavior for initial-boundary value problems for 
multi-term time-fractional diffusion equations with positive constant coefficients were studied by 
exploiting several important properties of multinomial Mittag-Leffler functions. m generalized 
the results in El by considering the counterpart with spatial-variable positive coefficients. 
In [12], the asymptotic behavior of the solutions to a more general time-fractional diffusion 
equations with distributed order derivatives were investigated. For other results for the time- 
fractional diffusion equations, we refer to m for a maximum principle, m, m for the initial¬ 
boundary value problems. 

In this paper, we are devoted to the long time asymptotic behavior for the space-time frac¬ 
tional diffusion equation. To the best knowledge of the authors, a long time asymptotic estimate 
for space-time fractional diffusion equation with nonpositive potential is not yet established. 

For later use, we introduce some notations. We first define the homogeneous Sobolev space 
H 1 (R d ) with 7 g R to be 

H^(R d ) := («g5': u(0) = 0, |£pu(£) G L 2 (R d )} 

equipped the norm 11u11 jj-y(R d ) := || || ^ 2 (]Rd) 7 and define the non-homogeneous Sobolev 

space FP(R d ) to be 

:= {u G S' : (1 + |£| 2 ) 3 w(£) G L 2 (R d )} 

with the norm |M| ff 7 (R<q := ||(1 + |£| 2 )^(£)|| L2(Rd) . 

We are ready to state our first main result. 

Theorem 1.1. Let 0 < a < 1 and 0 < 0 < 2 be given constants, a G 
p G -^“(iP) with p(x) <0 in R. d . Then 

I. ( t \n < f Ct ~* IMU 2 (R d )) rf </3 

~ \ct-“ ( ||a|| L1(R - ) + ||o|| ia(R - )) , d> 0 ’ 

is valid for the solution u to the Cauchy problem ED- 

Moreover, if the potential p satisfies p < —So, where So is a positive constant. Then there 
exists v(-,t) G #2 (R d ) be the solution of the elliptic problem 


LplP) n L 2 (R d ), 

t > 0 


(—A) 2 v(x, t ) — p(x)v(x, t ) 


tJt - Ta{x), x G R d , t> 0, 

F(l-a) 


(1.3) 
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which has the same asymptotic behavior as u, in the sense that 

IK-,0 -v{-,t)\\ H !3^ = 0{t- 2a )\\a\\ L 2 {Rd) , ast^-oo. (1.4) 

In the recent work [9] and m, by the use of Fourier analysis techniques, the asymptotic 
estimates of solutions to space-time fractional diffusion equations with potential p = 0 were 
studied. Their methods are based on the explicit representation formula of solutions, and thus 
cannot be used for the space-time fractional diffusion equations with non-zero potential. 

The above theorem shows that the decay rate of u is t~ a at best provided p < —<5o < 0. 
More precisely, we have the following statement: 

Corollary 1.2. Let 0 < a < 1 and 0 < /3 < 2 be given constants, a G L 1 ( R d ) H L 2 (R. d ), 
p G L°°(S . d ) with p(x) < —Sq in R d , where <5o is a positive constant. Moreover we suppose that 

l|u(-,f)|| p = o(t~ a ), as t — > oo. 

Then u(x , t) = 0 for all x G and t > 0. 

Now we turn to considering the asymptotic behavior of solutions to the space-time fractional 
diffusion equations in a bounded domain. Actually, we can consider a more general symmetric 
elliptic operator of second order such that A is positive, that is, there exists a constant So > 0 
satisfying ( Au,u) > < 5 oIM|| 2 (q) for all u G D(A). We consider the following initial-boundary 
value problem 

dfu = — A^u + p(x)u, x G 12, t > 0, 

< u(x,0) = a(x), x G 12, (1-5) 

_ u(x,t) = 0, x G <90, t > 0. 

We have 

Theorem 1.3. Let 0 < a < 1, 0</3<2 and 0 < 7 < ft be given constants, a G L 2 (Ll), 
p G T°°(0) with p{x) < 0 in O. Then 

\\u{-A)\\H-y(n) <Ct-° L \\a\\ L 2 {n) , f >0 ( 1 . 6 ) 

is valid for the solution u to the initial-boundary value problem (L3. 

Remark 1.4. For a G L 2 (0), as is proved in Section 3, in general we have 

u(t)<£H\ O) if 7 >/3. (1.7) 

Moreover the asymptotic behavior in time is determined only by the order of the fractional 
t-derivative. As Ob shows, the order of fractional x-derivative influences only the spatial 
regularity. 

The rest of this paper is organized in three sections. In Section 12.11 we investigate the 
analyticity of the solution to the Cauchy problem C3. By the Fourier-Laplace argument, 
we obtain an integral equation of the solution, and based on this, we further show that the 
solution can be analytically extended to the sector Sg 0 := {z G C \ {0}; |arg^| < d 0 } (f < 
do < min{ 7 r, ^}). In Section T2.21 we complete the proof of Theorem 1.1. In Section [3j we 
apply the eigenfunction expansion and the property of the Mittag-Lefflcr function to derive the 
asymptotic behavior of the solution to the initial-boundary value problem CD. Finally, the 
concluding remarks are given in Section [Tj 
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2 Asymptotic estimate in W l 

In this section, we mainly investigate the long-time asymptotic estimate of the solution to the 
Cauchy problem ed- Our proof is divided into two parts. In part 1, we will give an integral 
equation which is equivalent to ED- Based on this integral equation, we show the analyticity of 
the solution. In part 2, according to part 1 and the Laplace transform argument, the long-time 
asymptotic behavior is easily proved. 

2.1 Analyticity of the solution 

In this subsection, we are to give some properties of the solution u to the Cauchy problem ED, 
such as the continuously depending on the initial data, the analyticity, etc., which are useful 
for the proof of Theorem ll.il For this, we start from the following observations and notations. 
We first claim that the function 

w(x,t) = [ S a ,i(-|CI / 3 i“)wo(C)e 27 ri? ' a: d^ =: S(t)w 0 , x £ R d , t > 0 , ( 2 . 1 ) 

jR d 

where E at i(-) is the Mittag-Leffler function, see, e.g., [T9] , solves the Cauchy problem 

| d?w(x,t) = -(-A)%w(x,t), 0 <t<T, 

\ w(x, 0 ) = wq(x), x £ R rf . 

In fact, we can follow the Laplace-Fourier transform used in m to derive the above statement. 
By a direct calculation, we find 

S'(t)a(x) = - [ (2.3) 

J R d 

Based on the formula (12.11) of the solution to the Cauchy problem (12.21) . by regarding the term 
p(x)u(x , t) as a source term, we find an implicit form of the solution u to the Cauchy problem 
ED for t £ (0,T), say, the following integral equation of u, 

pt 

u(t) = S(t)a— / (— A)-%S'(t — r)(prt(r))dT, (2-4) 

Jo 

here we omit the spatial variable x, and we call a function u satisfies the above integral equation 
as the mild solution to the Cauchy problem ED- 

Lemma 2.1. Let 0<a<l,0</3<2 and 0 < 7 < /3 be given constants, a £ L 1 (R d ) DL 2 (M <1 ) 
and p £ L°°( R d ). There exists a unique mild solution u £ U 7 (R d ) to the Cauchy problem (11.11) . 
Moreover, the solution u : (0 ,T) —> H 7 (SL d ) is analytic and can be analytically extended to the 
sector Sg 0 , and the following estimate 


\H.,z)\\ HHRd) <Ce CT \z\ ^“|M| L 2 (K d), 0<7 </3, (2.5) 

is valid on a sector Sg 0t T := {z £ C \ {0} : | argz| < 0q, \z\ < T}, where ^ < 9 q < min{7r, 
and the constant C > 0 is independent of a, z, but may depend on d, a, /?, 7 , p, 9 q. 

Proof. After a change of variables in the integral equation (12.41) , we find 

C 1 

u(t) = S{t)a — t / (—A)~2 S' (t(l — T))(pu(t,T))dT. 

Jo 
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We now extend the variable t to the sector Sg 0 , and setting uq = 0, we define u n +i (n = 0, 1 , • • •) 
as follows: 

c i 

u n+1 (z) = S(z)a - z / S'(z(l - r))(pu n (zr))dr. (2.6) 

Jo 

We first claim that 

||Wi+l(-, z) — U n (-, 2)||z,2( R d) < C + i) kHMU 2 ^), z G Sg 0 . ( 2 . 7 ) 

In fact, first for n = 0, noting that the Mittag-Leffler function E atP (—z) (a > 0, p G R) is 

bounded on the sector Sg 0 , see, e.g., [IS] and 22 j we conclude from the Plancherel identity 
that the following estimates 

ll u l(‘> z ) — u o(-, ^)lli 2 (R‘ J ) = ll'S'(2:)a|| i 2( R d) = ||^a,l(—|CI^ 2 “)o(?)llL 2 (R rf ) < C|| a llL 2 (R d ) 

are valid on the sector Sg 0 . Now we assume that the statement (12.71) is valid for n — 1 G N, that 
is, 

rind ™ -1 

IM', Z) - U n - 1(-, ^)|| L 2 (Rd) < C _ + |~[(^- 1 )a|| a ||^ 2(Rd) , z G Sg 0 . 

For n G N, again that the Plancherel identity and the boundedness of the Mittag-Leffler function 

used in the evaluation for n = 0 yield 


r 1 

z / (-A ) -2 S'(z(l - T))p(u n - u n _i)(zr)dr 

J o 


||u n+ i(z) - u n (z 

= \z\ [ (~(! - 'r)) a_ 1 -Ea,c t (-|?| / 3 ( 2 (l - T )) a )(p{u n - Mn-l)(2T)) A dr 

Jo 

<C\z\ a [ (1 -t) q_1 ||(m„ -u n _i)( 2 r)|| L 2 (R d)dT. 

Jo 


L 2 ( R d ) 


L 2 { R d ) 


By the inductive assumption, we see that 


r(a ) n_1 

||un+i(^) - u„(2)||x,2(Kd) < - l)a + 1) 


r(a)r((n — l)o: + 1). 


r(na + 1) 


which proves the statement m - 

Now taking the operator (—A) 2 (0 < 7 < /3) on both sides of (12.61) . for n = 0, by noting 
the definition of u n and the Plancherel identity, we derive 

|| (—A) 2 (rti(z) — uq(z)) ||l 2 (r<j) = ||(—A) 2 S'(2:)a||i2( R£ i ) = 


\^S(z)a(0 


C 2 fR d 'l 


From the definition of the operator S(t) and the asymptotic expansion of the Mittag-Leffler 
function, see, e.g., Theorem 1.6 in m, it follows that 


||(—A)2( Ul (-) - M 0 (z))|| i2(Rd) < C 


igr 

i+i£H*i a 



Moreover, by noting that 


we have 


l + \Z\P\z\<* 


< C with 7 < /3, 


( 2 . 8 ) 


II(—A ) 2 ( 1 / 1 ( 2 ) — M 0 (z))||i 2 (R d) < 


\m ia 


i(OII- 


< C\z\ ? “||a||L 2 (R' 1 )- 

L 2 (R d ) 
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For n = 1,2, • • •, in view of the formula (12.31) and the asymptotic estimate of the Mittag-LefHer 
function, it follows that 


< C\z 


-A) ^ (u n+1 (z) - u n (z 

1 KPdsKl-T ))*- 1 




[ (-A) S'(z(l - r))p(u n - u n -i)(zr)dT 
Jo 

I (p(u n - u n -i)(zr )) A \dr 


L 2 (R d ) 


L 2 (R d ) 


Again the use of the inequality (12.81) and the assumption p £ L°°(R d ) imply 


(-A ) 2 (u n +i{z) - u n {z 


<C\z\ 




/V -r) c 

Jo 


(W)—1| 


- u„-i)(OT)|| L 2( R d)dr 


for z £ Sg 0 . Consequently, substituting m into above inequality, we can prove 


CY(a) n \z\ 


*(n+i-?) f i 


ll( —A) 2 (u„+i(z) - W n (2))|| i 2 (K d ) < - a + ^ 

Hence 


(i_ r) «(i-3)-i T - dT || a || i2(R£!) . 


(-A ) 2 (u n+ i(z) - u n (z 


C , r(a) n r(a(l - %)) , „ 


It is not difficult to see from the inductive argument that for each n £ N, u n (-,z) is analytic 
in sector Sg 0 under the norm I? 7 (R d ) with 0 < 7 < /3. Moreover, the classical asymptotics 


T{rf) = e v ri n 2 (2tt) 2 ^1 + O J J as 77 —* +00 
(e.g., Abramowitz and Stegun [T], p.257) imply that the series 

OO 

ll(-A)^ (u n+1 (-,z) - u n (-,z) 

n—0 

and 

OO 

ll U ™+l('^) - U n (-,z) 

n —0 

converge uniformly in any compact subset of Sg 0 , therefore 

OO 

u(-, z) := E(Vi(', z) - «„(•, z)) : Sg 0 ->• H' 1 (R d ) 


n—0 


is analytic. Finally, we see that u(-,t ) (i £ (0,T]) is just the solution to the Cauchy problem 

(EH). 

From the integral equation (12.41) of u, similar to the above argument, we see that 

r\z I 

ll'w(-^)IU 2 (R d ) < C'l|a|U 2 (Rd) +C / (| 2 | - r) a_ 1 ||it(-,z)|| L 2 (R d ) dr. 

Jo 

Thus by the generalized Gronwall inequality (see e.g. Lemma 7.1.1 in [ 6 ]), we find 

\\ u {-: ~)IU 2 (R d ) < Ce CJ ||a||i 2 ( R d), 0 £ Sg 0> T, 

combining with the integral equation (12.41) . hence that 

||(—A ) 2 u(-, z)|| L 2 ( R d) < C'e c " r |^| ?a ||a||i 2 ( R d), z £ Sg 0i T- 
This completes the proof of Lemma 12.11 □ 
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2.2 Proof of Theorem 11.11 

In this part, we are to finish the proof for the asymptotic behavior of the solution u to the 
Cauchy problem (HUD- We first investigate the property of the Laplace transform £u of the 
solution u. For this, we apply the Laplace transform 

pQG 

(Cf)(s) := / f(t)e~ st dt 
Jo 

on both sides of the differential equation in (ED and use the formula 

C(d?f)(s) = s a (Cf)(s)-s a - 1 f( 0+), 
to derive the transformed algebraic equation 

s a £u(-s) =—(—A) 1 ^£u(-; s) + p(x)£u(-; s) + s a ~ 1 a, s £ C, Res>M, 
where M > 0 is a sufficiently large constant. 

We check at once that £u :{s£C: Res > M} —> (R d ) is analytic, which is clear from 
the property of the Laplace transform. Moreover, we claim that £u(x\ s ) (Re s > M ) can be 
analytically extended to the sector {s £ C \ {0} : j args| < 7 r}, that is, 

Lemma 2.2. Under the assumptions in Lemma \2.1l the Laplace transform £u of the unique 
solution u to the Cauchy problem ED can be analytically extended to the sector {s £ C \ {0} : 

I arg s| < 7 r}. 

Proof. By Lemma 12.11 it follows that u(z ) is holomorphic on the sector |argz| < 0 < ^ 
under the norm in H 1 ( R d ) (0 < 7 < /3), such that (12.51) . we define /(s) = £u(s),s > 0. 
For cj) < 9, consider the path 7 #^ = re 1 ! 0 ’^, [r, i?], Re 1 ^’^, [R, rje 1 ^; by the Cauchy theorem 
f u(z)e~ sz dz = 0. i.e. with r —\ 0, R —>• 00 we obtain 

J'YR.r V ' ’ ’ 

f-oo nOO 

f(s) = / u(t)e~ st dt = / u(te^)e~ ste ^e'^dt, s > 0 . 

Jo Jo 

Since </> £ (—9,9) can be arbitrarily chosen, it follows that f(s) can be analytically extended 
to the sector {s £ C \ {0} : | args| < 9 + (V0 £ (0, f)). This complete the proof of Lemma 

O “ □ 


See e.g., Theorem 0.1, p.5 in [23] for a similar argument. Moreover, we have the following 
lemma about the evaluation of the Laplace transform £u on the sector Sg 0 . 

Lemma 2.3. Let 0 < a < 1 and 0 < /3 < 2 be given constants, a £ L 1 (IR d ) n L 2 (M. d ), 
p £ L^R^) with p(x) < 0 in R d . There exists a constant C depending on d, 9 q, a, [J such that 


l|£u(-;s)|| • | 


< 


C|s|= 1 ||«||z, 2 ( R d ) , d</3 

C\s\ a 1 (||a|| i i( R d) + ||a||i 2 ( K d)), d > ft 


Proof. From the above argument, we see that 


Vs £ Sg 0 . 


(2.9) 


(—A) 2 £ u (-; s) + (s a — p(x))£u (-; s) = s a 1 a, s £ Sg 0 . (2-10) 

Firstly for any fixed s = re' e £ Sg 0 , we define a bilinear operator B[«L>, \R; s] : H%( R d ) x 
Rv(R d ) by 


B[<h,T;s] 



((-A)^(x)) dt(x) + (s a 


p(a;))$(a;)'I , (a:)da;, 
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where '3/ denotes conjugate of 4'. Consequently, for any s £ Sg 0 , (12.101) can be written as 
B[£u,*;s] = VTei^(R d ), s € S 9o . 

Furthermore, taking 'F = Cu in the above equality, from the Parseval identity and Fourier 
inversion transform, we have 

B[£rt, Cu\ s] = / ((—A) 2 £u(£;s)\ s)d£ + / (s“ — p(x))Cu(x] s)Cu(x\ s)di 

Jm. d ' ' J«. d 

= [ (|£|^£u(£;s)) (|^|^£u(^;s)) d^+ [ (s a - p(x))\Cu{x; s)\ 2 dx 
JWL d V J V 7 JR d 

= f (— A) 1 *£u(x; s)(—A)^Cu(x- s)d:r + f (s a — p(x))\Cu(x; s)| 2 dcc. 

Jm. d J R d 

In view of cos (ad) > 0 for 0 < 9 < and the assumption p < 0, it follows that 


ReB[£u, Cu\ s] > (—A) 4 £u(-; s) 
In the case of d > /?, we conclude from (12.111) that 


L 2 (WL d ) 


■q«n| J Cu(-;«)|| ia(R - ) . (2.11) 


||£w(- ;s)||^| < C\Q[Cu, Cu\ s}\ < Cjs|“ |(a, Cu (-; s)) i 2 (R d ) |. 


We note that 


(a,Cu{- ;s)) L 2 {Kd) = J ^a{x)£u{x;s)dx = J ^ (|£| 2 a(0) (|£| = £u(£; s)) d£, 

from the Holder inequality, we see that 




'R d 


|£|-*|a(0| 2 d£] / |£|* Cufcs) 




d? , 


which implies 


ii J c«c- ; a )ii ^ <cr|«r 


Ier^|a(0| 2 d^ 


s £ Sg 0 . 


Breaking the integral on the right-hand side of the above estimate into {|£| < 1} and {|£| > 1}, 
in view of ||a(||x,° o (m d ) A IHIl^r 1 *), it follows that 

/ kl^|S(0| 2 de < \\a\\ 2 LHMd) f lei^dC + [ |a(£)| 2 d£ < C(\\a\\ 2 Lim + ||a||£ 2(Rll) ), 

jR d V ’J Iflcl J\e\>l V ’ 


hence that 


l|£ w (' ; S )ll^§( R d) — C| s r 1 (ll a llL 1 (K d ) + IHlL 2 (R d ))> s£Sg 0 ,d>p. 


Next for d < /3, noting that |H|z, 2 (M ,q||£u|| L 2 (Rd) < ^|s| 1 ||a||| /2(R d ) + e 0 |s|||£u||| 2(Rd) with 
sufficiently small constant eo > 0, we see from (12.111) that 


l|Au(- ; s)||^| (Rd ) <C|s|= 1 1|« 


S £ Sg 0 , d < f3. 


The proof of Lemma T2.31 is completed. 


□ 




















Now let us turn to finishing the proof of the long-time asymptotic behavior of the solution 
u to the Cauchy problem (dis¬ 
proof of Theorem 11.11 By the Fourier-Mellin formula (e.g., [23]), we have 

^ pM-\-ioo 

u{x,t) = ■— / Cu(x; s)e st ds. 

JM—ioo 

From Lemma we see that the Laplace transform £u(x; s) of the solution to the equation 
(12.10H is analytic in the sector Sg 0 . Therefore by the Residue Theorem (e.g., [21] ). for t > 0 we 
see that the inverse Laplace transform of Cu can be represented by an integral on the contour 
7(e, 0o) defined as {s £ C : args = 9q, |s| > e} U {s € C : |args| < 9 0l |s| = e}, that is 

u[x,t) = / £u(x; s)e st ds, 

2lT1 Jj(e,9 0 ) 

where the shift in the line of integration is justified by the estimate m • Moreover, again from 
the estimate ED , we can let e tend to 0, then we have 

u(x,t) = —: / £u(x; s)e st ds. 

2lrl J 7 (0,ff 0 ) 


Thus 


||u(-,i)||.£ <c/ ||£u(-; s) II . 3 |e st ds|, (2-12) 

11 v - J j(0M 11 Hz (R d ) |! v ; 

hence combining (12.91) with (12.121) gives 

/*oo 

\\u{-,t)\\ <C / r^ _I ||a|| L 2 (R d ) e rtcos9o dr < |M| L 2 (R d), t > 0, d < /3, 

H 2 ) J 0 


and 


IK- ’ < )ll i j# (R d) ^ C J r “ 1 H a lli 2 (R- i ) er ' tCOSeodr < ct “(ll«llz, 1 (K‘') + l|a||L 2 (Rd)), t > 0, d > (3. 

Now let us turn to consider the asymptotic behavior of the solution u in the case of p < —do, 
where do is a positive constant. We repeat the above argument to derive that the Laplace 
transform Cv of the solution v to the equation CD such that 

(-A)^Cv(-;s)—p(x)Cv(--,s) = s a ~ 1 a, s £ Sg 0 . (2-13) 

We define a bilinear operator Bi [<I>, \R] : H^(R d ) x H^(R d ) —> C corresponding the equation 
(12.131) by 

Bi[$,’F]:= J ^(— A)^$(x)j *F(x) — p(x)<h(x)’I'(x)dx. 

By noting p(x) < —do < 0 in similar to the proof used in Lemma T2.31 we find 

||£u(-;s)||^| (Rd) < ReBi[£u, Cv\ < C|s| a_1 |ja|| L 2 (R d) ||£u(-; s)|| L 2 (R d), s £ Sg 0 . 

Noting that Cu — Cv satisfies the following problem 
(—A )%(£u(-;s) - Cv(-;s )) + (s“ - p(x))(£u (-; s) - £v(- ;s)) = s a £v(-;s), s £ Sg 0 . 
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Then again by the argument used in the evaluation of (12.91) . we deduce that 
||£u(-;s) -£«(•; s)|| H | (Rd) < C^WCv^ ; s)|| L 2 (R d). 


Consequently 


\\u{- ,t) -v(-,t)\\ H z 


< C 


'7(0, So) 


l|£(«(- ;s) ^ *>(•;«) 




*ds| 


pOO 

<c r 2 “ _1 ||a|| i 2 (M <i ) e rtcoseo dr < Ct _2a ||a|| L 2 (Rd) , t > 0. 

Jo 

This completes the proof of Theorem 11.11 □ 

Proof of Corollary 11.21 The relation 

\\t a u(-,t) - = t a O(t 2 “)|H|i,2(m<*) -> 0, t -» 00 

follows easily from the estimate m in Theorem 11.11 By the assumption of Corollary 11.21 
||i“u(-,t)||^|^ d ^ —» 0 as t —> 00 , which implies the relation t a v(-,t) —> 0 as t —> 00 under the 

norm || • ||^£^ d ^. By noting that v solves the equation m, it follows that 


Bi[v(-,t),ip\ = r _ (a, (p) L *(R*), p€H* (IT), 

which implies 

Bi[tv(-,t),ip\ = (a,<P)L 2 (R*), V e H%( R d ). 

Letting t —>• 00 , we have 

= 0 in R d . 


T(l-a) 

In this case, the Cauchy problem m has only the trivial solution, e.g., u(-,t) = 0 in and 
t > 0. This completes the proof of Corollary 1 1.2 1 □ 

Remark 2.4. The existence of the unique solution v to the equation (O can be easily shown 
by using the Lax-Milgram theorem. Indeed, we define a bilinear operator Bi[$,4'] : H 2 (R d ) x 
H* (M d ) —>• C corresponding the equation <ra by 


B 


[$,4']:= f ((— A) 2 $(x) S ) 'L(x) — p(x)$(x)4'(x)da 


Similar to the proof used in Lemma. 1 2. ,91 we see that 


|Bi[$,®]| <C||$|| e Ill'll e 

1 L ’ Jl - 11 (R^) (R<J)’ 

moreover, by noting that p(x) < —<5o < 0 in R d , we have 

Then the Lax-Milgram theorem shows that there is a unique solution v G H% (R d ) to the equation 
(ED which satisfies 


Bi[v(-,£), ip] 


t~ a 

P(1 _ a ) ( a ’^ 2 (R d )> 


p G H%( R rf ). 
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Thus by taking (p = v, we have 


Ct~ a 

IK’011% - C' B ib(-,t) J w(-,t)] < —- -\(a,v{-,t)) L 2 {Rd) \. 

H 2 (Rd) 1 (1 — a) v ’ 

Hence 

lk(-,i)|| ff 4 (Rd) < Ct- a \\a\\ L 2 {Kd) . 

3 Asymptotic estimate in a bounded domain 

In this section, wc will prove Theorem 1.3. 

We first show two useful lemmas. 

Lemma 3.1. The operator A% — p is positive and self-adjoint. 

Proof. Since A% is positive and self-adjoint, noting that the perturbation by p is bounded, we 
see by Theorem 4.3 (p.287) Kato [ 8 ] that A% — p is self-adjoint. The positivity follows from 

p < 0. □ 

Lemma 3.2. There exist constants C±,C 2 >0 such that 

Ci\\Ai a||i 2 (n) < \\ — A? a + pa\\ L 2 (Q,) < C 2 \\A^a\\ L 2 ^ 

for all a £ D(A%). 

Proof. First we note that ||a||z. 2 (£ 7 ) < C\\A %for a £ D(A^). Therefore 

|| — A? a +pa\\ L 2 ^ < \\A^ a|| L2 (Q) + C'||a|| i 2 (Q) < C||^4' 2 'a||L 2 (n) 

for a £ D(A^), which is the second inequality. Next since A% — p is positive, 0 £ p(—A% +p): 
the resolvent. Therefore (—A? +p) _1 : L 2 (fl) —> L 2 (fl) is bounded. We set b = (— A 2 + p)a, 
and a = (—A% +p)~ 1 b. Then 

P^a|| L 2 (n) = \\A^(-A^ +p)~ 1 b\\ L 2 {n) = || (A% - p)(-A % +p)~ 1 b + p{-A% +p)~ 1 b || L 2 (n) 
<IHU 2 (n) + +p) 1 ^IU 2 (n) < IHU 2 (n) + C'l|&IU 2 (fi), 

that is, 

ll^allz^n) < (1 + C)\\(-Ai +p)a|| i 2( n ), 

which is the first inequality. Thus the proof of the lemma is completed. □ 

In view of Lemma 1X31 we apply the Heinz-Kato inequality (e.g., Tanabe (31]) and have 

D{{-A% +p) e ) = D(A% e ) 

and 

C'i||A 2 ' 6l a|| i 2 (Q) < ||(— A~ +p) e o|| i 2 ( 0 ) < C 2 \\A^ e a\\ L 2 ^, a £ D(A^ 8 ) 1 0 < 9 < —, 
that is, 

Ci\\Ai a\\ L 2(Q) < ||(—d 1 +p) ? a||L 2 (n) < C 2 \\Aia\\ L 2^) (3.1) 
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and 


D((-A% +p)%) = D(A%), 0 < 7 < 2. 

Let X n be the eigenvalue of A% — p: 0 < Ai < A 2 < A 3 < • • •, where A j appears the same times 
as the multiplicity. Let ip n , n 6 N, be an orthonormal eigenvector for X n . Then {</3 n }rieN is an 
orthonormal basis in L 2 (f2). By (13.11) . we have 


Ci|l a llffT(n) < 


X n (a, Pnj 1 


< C 21 |«|| 


(3.2) 


for a € D(Al ) and 0 < 7 < 2. 

Now we proceed to 

Proof of Theorem I1.3L Similarly to Sakamoto and Yamamoto [22], we represent the solution 
u by 

OO 

u(t) = y^(a, Pn)E atl (-X n t a )ip n 

n— 1 

in C((0,T\;D(-A% +p)). By (|3l2l) . we estimate ||u(f)||jj- 7 (n): 


\\um 2 HH n) <cY. X ^(a,p n ) 2 \E aA (-X n n^ 


<C^(a ,^) 2 


n= 1 


n =1 
' X 
Xn 

1 + x n t° 


= ct p 


n— 1 


( (A n t a )^ 

Y1 + A n t a 


2 


Let 7 < /3. Then j— is monotone decreasing for p > -gbz- For sufficiently large t > 0, we 

can assume that Aii“ > -p— ■ Therefore 

p -7 

(A „t“)* A ft? „ a, a 

\ n , < - L _- < Ctlrt - a 


1 + A n t a ~ 1 + Ai t° 


that is, 


< Ct " a YXa, pn) 2 


n —1 


for large t > 0. For 7 = /?, since < 1, we can directly see the conclusion of the theorem. 

Thus the proof of the theorem is completed. □ 

Proof of (11.711 . Let 7 > /?. By Theorem 1.4 (pp.33-34) in Podlubny [T9] . we have 


E a ,i(—X n t a ) = 


+ 0 


X n t a r(l — a) \A 2 f 2 


We set j = 1 + 5 with small S > 0 and set itjv(f) = J2n- i(°; fn)E a ^(—X n t a )p n . Then 

N 

il^W 111/7(0) = X2 n +2 \ a ^n) 2 \E aA (-X n t a )\ 2 


n =1 


N 


= A 2 + 25 (a, ip n ) 2 


1 1 


+ 0 


A„t“ T(1 — a) \A 2 f 2 


1 
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N 


= 5>i;«W) 2 


n=1 
N 


1 1 

A n t a T(1 - a) 


o 


f 1 


\K t3c 


Y^x 2 n \a,ip n ) 2 ^jt 2 “(f(T^^y) ~ c Y^ x n x ( a ^n) 


> 


2 t~ 3a . 


Here we note 


N 


E \ 2 S- 1 / \2.-3a 

"n K a j l Pn) t 


< £ A“-W) V>« 

n=l 


< OO 


because S > 0 is small so that 2(5 — 1 < 0. Thus we verified that a §£ H 2S (Q) implies 
hnpv _>. 00 Y^n= i( a > Vn) 2 )^ = oo and so u(t) ^ □ 


4 Conclusions and remarks 

In Lemma 12.11 we proved that the solution u(-,t) : (0, T) —»• iL 7 (R d ) (0 < 7 < /3) is analytic 
in time t £ (0 ,T) and that the short-time estimate ||u(-,t)||# T ( R d) < Ce CT t~P a \\a\\ L 2 ^ R d' ) as 
t — > 0. Moreover, u : (0 ,T) — > H 7 (R d ) can be analytically extended to the sector Sg 0 = {z £ 
C \ {0} : |argz| < 0 O } with 9 0 such that ^ < 9 0 < min{ 7 r, j-}. However, we do not know 
whether the solution can achieve the sharp spatial regularity H^(JH d ). 

In Theorem 11.11 we discuss the long-time asymptotic behavior of the solution to the Cauchy 
problem (O). We show that the decay rate of the solution is dominated by t 2 in the case of 
d < f3 and by t~ a for d > /3 under the assumption that the potential p is in the space L°°(R d ) 
and such that p < 0 in M. d . However, from m and the property of the Mittag-Leffler function, 
it is not very difficult to see that the solution u to the Cauchy problem O) admits the sharp 
estimate 

]K-,i)||^f (R(i) < Ct~% IMIz^d), d > 1, 0</3<2, as t 00 

and 

||w('4)ll^ (Rd) < ^““(HoHiifRd) + || a llL 2 (R d ))) d > /3 , as t —> 00 

provided that p = 0 in that is the estimate in Theorem 11.11 However, we do not intend 
to discuss whether the first result in Theorem 1 1.1 1 is sharp in the case of p < 0. If we assume 
p < —5o < 0 , the decay rate is shown to be t~ a and sharp. 

In Theorem ll.31 the long-time asymptotic behavior of the solution to the initial-boundary value 
problem is investigated. By the property of the Mittag-Leffler function and the eigenfunction 
expansion, the solution to (1.5) decays as t~ a which is different from the case of unbounded 
domain, where in the case of H = R d the decay rate is shown to be dependent of the relation 
between /3 and d. 


Acknowledgement 

The research of Xing Cheng is supported in part by “the Fundamental Research Funds for 
the Central Universities” (No.2014B14214). Zhiyuan Li thanks the Leading Graduate Course 
for Frontiers of Mathematical Sciences and Physics (FMSP, The University of Tokyo). The 
works was partly supported by A3 Foresight Programs “Modeling and Computation of Applied 
Inverse Problems” by Japan Society of the Promotion of Science. 


13 





References 


[1] M. Abramowitz and I. A. Stegun, Handbook of Mathematical Functions with Formulas, 
Graphs, and Mathematical Tables, Dover, New York, 1972. 

[2] Z. Chen, M. M. Meerschaert and E. Nane, Space-time fractional diffusion on bounded 
domains, Journal of Mathematical Analysis and Applications 393 (2012), 479-488. 

[3] A. Friedman, Partial Differential Equations of Parabolic Type , Dover Publications, 2008. 

[4] R. Gorenflo and F. Mainardi, Random walk models for space-fractional diffusion processes , 
Fract. Calc. Appl. Anal 1 (1998), 167-191. 

[5] R. Gorenflo, F. Mainardi, D. Moretti, G. Pagnini and P. Paradisi, Discrete random walk 
models for space-time fractional diffusion , Chemical Physics. 284 (2002), 521-541. 

[ 6 ] D. Henry, Geometric Theory of Semilinear Parabolic Equations, Springer-Verlag, Berlin, 

1981. 

[7] H. Jiang, F. Liu, I. Turner and K. Burrage, Analytical solutions for the multi-term time- 
space Caputo-Riesz fractional advection-diffusion equations on a finite domain, Journal of 
Mathematical Analysis and Applications, 389(2012), 1117-1127. 

[ 8 ] T. Kato, Perturbation Theory for Linear Operators, Springer-Verlag, Berlin, 1976. 

[9] J. Kemppainen, J. Siljander and R. Zacher, Representation of solutions and large-time 
behavior for fully nonlocal diffusion equations, arXiv:1505.02803 

[10] K.-H. Kim and S. Lim, Asymptotic behaviors of fundamental solution and its derivatives 
related to space-time fractional differential equations, arXiv:1504.07386. 

[11] Z. Li, Y. Liu and M. Yamamoto, Initial-boundary value problem for multi-term time- 
fractional diffusion equation with positive constants coefficients, Appl. Math. Comput. 257 
(2015), 381-397. 

[12] Z. Li, Y. Luchko and M. Yamamoto, Asymptotic estimates of solutions to initial-boundary- 
value problems for distributed order time-fractional diffusion equations , Fract. Calc. Appl. 
Anal. 17 (2014), 1114-1136. 

[13] Z. Li and M. Yamamoto, Initial-boundary value problems for linear diffusion equation with 
multiple time-fractional derivatives, arXiv:1306.2778 

[14] Y. Luchko, Maximum principle for the generalized time-fractional diffusion equation, J. 
Math. Anal. Appl. 351 (2009), 218-223. 

[15] Y. Luchko, Initial-boundary-value problems for the generalized multi-term time-fractional 
diffusion equation, J. Math. Anal. Appl. 374 (2011), 538-548. 

[16] F. Mainardi, Y. Luchko and G. Pagnini, The fundamental solution of the space-time frac¬ 
tional diffusion equation, Fract. Calc. Appl. Anal. 2 (2001), 153-192. 

[17] F. Mainardi, P. Paradisi and R. Gorenflo, Probability distributions generated by fractional 
diffusion equations, arXiv:0704.0320. 

[18] A. Pazy, Semigroups of Linear Operators and Applications to Partial Differential Equa¬ 
tions, Springer-Verlag, Berlin, 1992. 


14 


[19] I. Podlubny, Fractional Differential Equations, Academic Press, San Diego, 1999. 

[20] J. Priiss, Evolutionary Integral Equations and Applications , Springer Science & Business 
Media, 1993. 

[21] W. Rudin, Real and Complex Analysis (3rd), New York: McGraw-Hill Inc, 1986. 

[22] K. Sakamoto and M. Yamamoto, Initial value/boundary value problems for fractional 
diffusion-wave equations and applications to some inverse problems, J. Math. Anal. Appl. 
382 (2011), 426-447. 

[23] J. L. Schiff, The Laplace Transform Theory and Applications, Springer, 1991. 

[24] H. Tanabe, Equations of Evolution, Pitman, London, 1979. 

[25] M. Yamamoto, Carleman estimates for parabolic equations and applications, Inverse Prob¬ 
lems, 2009, 25(12), 123013. 


15 


